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A combined scattering mechanism involving the states of free electron-hole pairs (exciton contin- 
uum) and discrete excitons as intermediate states in the multi-phonon Raman scattering leads to 
(1) a strong increase of the scattering efficiency in the presence of a high magnetic field and to (2) 
an outgoing excitonic resonance: the two features are not compatible when only free pairs (leading 
to a strong increase of the scattering efficiency under the applied magnetic field) or discrete excitons 
(resulting in the outgoing resonance at the excitonic gap) are taken into account. 
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I. INTRODUCTION 

In a recent publication^ we have shown that the strong 
outgoing resonance observed in high order multi-phonon 
resonant Raman scattering (MPRRS) from polar semi- 
conductors can be explained when high energy interme- 
diate electronic states belong to the excitonic continuum 
(approximated by free electron-hole pairs, EHP) and only 
couple to the bound excitonic state at the last stage of a 
scattering process. The high probability of decay into the 
continuum strongly opposes the MPRRS mechanism in- 
volving discrete excitons as the only intermediate states 
for explanation of the observed outgoing resonance at 
the ground excitonic transition (see Ref. [j] and references 
therein). Cooled by the emission of a sufficiently large 
number of LO-phonons, the EHP binds into an exciton 
whose energy is not enough for LO-phonon-assisted de- 
cay. 

In this work we analyze the effects of a high magnetic 
field on the outgoing excitonic resonance considering, as 
in Ref. ^, the monomolecular creation of a cold exciton 
by the light-generated free EHPs which lose energy but 
preserve their spatial correlation through the interaction 
with LO-phonons. 
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FIG. 1. The diagram involving one discrete exciton in- 
termediate state contributing to the MPRRS efficiency in 
the range of outgoing resonance. Hollow circles represent 
photon-electron-hole pair interaction, bold circles correspond 
to the electron-LO-phonon interaction while the square ver- 
tices are shown for discrete-continuum transitions. Solid lines 
above (below) the dash-dotted line represent the electrons 
(holes) and horizontal lines stay for bound excitons. Bold 
dashed lines connected left and right hand sides of the dia- 
grams correspond to LO-phonons while the dotted lines rep- 
resent incident (on the left and right sides) and scattered (in 
the center) photons. 



II. MODEL 

The main contribution to the A^-th order MPRRS ef- 
ficiency follows from processes with one (Fig. and two 
(Fig. ||) bound excitonic intermediate states at the last 
stage of the elementary scattering process. Only these 
contributions correspond to the cascade of transitions 
where the bound exciton cannot decay into the EHP con- 



tinuum through the emission of LO-phonons. We assume 
to/j > > m e so that the hole energy is less than the energy 
of one LO-phonon and all phonons emitted by the EHP 
before its binding into a discrete exciton are emitted by 
the electron. 

We use the Landau gauge A = A(0, xH, 0) for a mag- 
netic field directed along the z-axis and the corresponding 
wave functions of free EHPs. Only the ground state of 
the bound exciton is taken into account for discrete exci- 
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to be used below. For £(s) = exp (— |s|) one finds £ = 1, 
77(a) = 1/(1 + a 2 14) and 6(a) = 2/(1 + a 2 ). The wave 
function of Eq. (|l) reduces to that of Ref. || when the 
Landau gauge is changed for the symmetric one. 



III. SCATTERING EFFICIENCY 



FIG. 2. The diagram with two discrete exciton intermedi- 
ate states contributing to the MPRRS efficiency in the range 
of outgoing resonance. The square vertices are shown for 
transitions between two states of discrete exciton and for dis- 
crete-continuum transitions. 



tonic intermediate states in the last stage of the process. 
According to Refs. ^ and [|, the bound exciton wave func- 
tion in a high magnetic field (a >> an, where a is Bohr 
radius in a zero magnetic field and an is the magnetic 
length, an = yhcfeH) can be written as 
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exp 
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x exp{i [(K x - (y/a 2 H ))R, 
+ $(r ± ,-K ± )+C(K x ,K y )}} 



KyRy 



(1) 



and 



$(r ± , -K x ) = (xy/a 2 H - K x x ~ K y y) (m e ~ m h )/2M , 
C(K x ,K y ) = a 2 H K x K y (m e - m h )/2M , (2) 

77i e (t«/i) is the electron (hole) effective mass, M = m e + 
TO/j, rj_(Kj_) = (a 2 H /H) [H x Kj_] and R, r are the center 
of mass and relative motion coordinates of an electron 
and hole. The longitudinal part ^\\k, of the exciton wave 
function can be written as 
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where £(s) describes the relative motion of the electron 
and hole along the magnetic field direction and satisfies 
the equation £(s = 0) = 1. The constant £ is determined 
by the normalization condition 



ds £ 2 (s) . 



We do not specify the exact form of £(s) (see Ref. ||) and 
introduce the two functions 
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The scattering efficiency can be written a: 
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e* sa e sf 3ei 7 e^ x S a jf3x{uJi,uJ s , m, n s ) , (4) 



where S a yp\ is the light scattering tensor of rank four, c 
the light velocity in vacuum, ni (n s ), e/ (e s ), k; (k s ), ui 
(u s ) are the refractive index, polarization vector, wave 
vector, and group velocity of the incident (scattered) 
light, respectively. Using diagrammatic techniques, sim- 
ilar to those of Refs. [TJ, || and|5[ we find for the contribu- 
tions of the diagrams in Fig. nT (Fig. 0) 
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where 00 = (e 2 /moC 2 ) 2 and jexcniO) is the inverse life 
time (broadening) of the exciton at the ground state 
with energy Em = hcum- According to the assump- 
tion mh » 777 e , the energy and the broadening of the 
hole have been neglected in all energy denominators. The 
quantity Ljv a m has dimensions of length and, for the di- 
agram in Fig, 
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where the index (3 designates the sequence of transi- 
tions made by the electron through Landau bands emit- 
ting successively N — 1 phonons. It represents the 
set of 77o, 77i, . . . , 77tv_i Landau numbers and indices 
7i,Z2, . . . ,ijv_i- Each index i may be or 1: it is zero 
when the electron does not change the direction of mo- 
tion along the magnetic field after the phonon emission 
and one when the sign of the velocity is opposite in the 
states before and after the phonon emission. In Eq. (0), 
Dm p represents the integral 
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where 
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and 



X i (K,K',x)= [x + a 2 H (K T K') 2 /2\ \ (10) 

lie), where 
LQ), the — (+) sign corresponds to 



fudgii = Eg + heH / (2uc) , where \i — m e mh/M and E g 
is the gap. In Eq. 
i = (i= 1). 

Note that ijy = when the direction of motion along 
magnetic field after emission oi N — \ phonons coincides 
with an initial direction. In this case, s = ^Z^Lj 1 i n 
is an even number. For odd s the direction of mo- 
tion is opposite and ijv = 1. The symbol (...) corre- 
sponds to the average over the directions of wavevec- 
tors qiX) Q2_L> • • • > QiV— l-L m the xy-pl&ne, when Xj — 



a 2 H q 2 ± /2 and q N ± = - Y^iL^ 1t-L • We used also 

AT-l 

Fjv^ = JJ (2-Yj/otULo) , 7j = 7e(n>j,Kj) , (11) 

where a is the Frohlich coupling constant. When 
7e(n, |fc z |) is determined by the interaction with LO- 
phonons, we find 



(2l\k' z \)- 1 dxB nn ,(x) 



x [x°(|A,|,|Ai|,aj)+x 1 (|**|,|Ai|,aj)] , 

= V 7 ^ + 2m e [w eff (n - n') - w LO ] /fi . (12) 

The sum over n' is limited by the condition that k' z has to 
be real. When all 70, 71, ... , 7jv-i are determined by the 
probability to emit an LO-phonon in a real transition, 
the substitution of Eq. (|l2| ) in Eq. ( 11 ) and multiplying 
the result by l N leads to Y/v/3 defined in Eq. (134) of 
Ref. |[ However, close to the resonance ui s = u>m, the 
electron emitted N — 1 LO-phonons occupies the state 
with the energy less than the energy of an LO-phonon. 
In this case jn-i is determined by some other weaker 
scattering mechanism: 
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This leads to the result: 
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The length A^p is defined as 



= fffp(z = 0) , fNp{z)=f ili >" i »- 1 {z). (15) 
For example, in the case N = 4 

f W (z) = [f ++ -(z) + f- ++ (z)]/2 , 
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J «0,"w-1 >4 [^ ,o9 ( a ||-^o) + <W_i,o© (fl||^iv-i) 



Aj = / m e 7j and 7^ = 7 e («.,■ , Aj ) . 

In ^-representation / +H (z) can be written as 

I roc 

/++ (z) = — / dpexp(ipz) 
[{I + i\ oP ){\ + i\ lP ){l ~ i\ 2P ){l - iX^p)]- 1 . (18) 
Finally, we have used the definition 
1 

c 3 

- 2S no , 5 nN _ 1 , e(a l] K )Q(a l] K N - 1 )}, (19) 

where 0(a) has been defined after Eq. (||). 

We proceed to calculate the contribution of the di- 
agram in Fig. ||. Since one of the intermediate states 
for the process of Fig. || corresponds to an exciton with 
K / 0, we need to comment on some details of the 
ground exciton dispersion E exc (K±,\K z \). The energy 
E exc (K±, \K Z \) can be written as 

E exc (K ± , \K Z \) = E gH + E(K ± ) + ti 2 K 2 /(2M) . (20) 

The function E(Kj_) in some limits can be found in 
Rcf. ||. For our purposes it suffices to note that E(K± = 
0) = — AEih is the exciton binding energy in a high mag- 
netic field and E\h = E g n ~ AEih- The contribution of 
the diagram in Fig. ^ is given by Eq. (|5|), where 



1 M I 

= au>LO — / axexp(—x) 

L Nb m, 
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K za = y/2M [lji - Lu gH ~ E{x)/h -{N- 1)lo lo ] ]% is 
the absolute value of the z-component of an exciton wave 
vector in the N-th real intermediate state, x = a 2 H K\/2, 
and T](a) has been defined after Eq. (||). The K±_ max 
is a maximum value of K± allowed by energy conser- 
vation, i.e., under the condition that K z q is real. At 
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variation of Kj_ from zero to, infinity the value E(x) 
changes from —AEw to zero.Q Therefore, in the range 
uji > ujgH + (N — l)u> XjOi w e have x max — > oo, whereas 
for luih + (N— 1)lu lo < < Lu gH + (N— 1)uj lo the val- 
ues of K± max and x max are determined by the equation 
fiLJi - E gH - (N - \)huj LO = E(K± max ). We used also 
the following definitions: 



Rn-i/3(K±,K z0 
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(22) 



P (K ± , K z0 ) = [T fi (K x ,K z0 ) + T (K±, -K z0 )] /2 
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indexes no, n\, . . 
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set 



(24) 
of 



The variables of integration in Eq. (|22|) are Xj = a| f g|_ L /2 

with a constraint qjv-i_i_ = — Xii=i 2 1i-L — The sym- 
bol (• • •} denotes the average over angles which determine 
the direction of vectors qi±, q 2 _i_, ■ ■ ■ , qAr-2±, K_l in the 
xy-plane. 



IV. APPLICABILITY LIMITS 

Let us discuss the applicability limits of the expressions 
for contributions of diagrams in Fig. [j] and Fig. |[ We 
assume that AEh < Tiojlo and consider four intervals 
for the laser frequency: 

(a) UiH + (N — l)bJLO < tOl < UJgH + (N - 1)uolo , 

(b) ujgH + (N - 1)uj lo < loi < loih + Nlo lo , 

(c) OJXH + NuJ LO < LOl < LOgH + Nuj LO , 

(d) uji > uj gH + Nuj lo ■ (25) 

The width of the intervals (a) and (c) is AEh/U and 
the one of interval (b) is [lolo — AEh/Ti). The outgoing 



excitonic resonance coincides with the border of intervals 

(b) and (c). 

Equations ( |) and (^J) for the contribution of the dia- 
gram in Fig. |l are valid when K^-i is real for n^-x = 
(see Eq. (ft)). This is satisfied within the intervals (b), 

(c) and (d). In intervals (b) and (c), the value of 
Kn{tin = 0) is pure imaginary and it is real in inter- 
val (d). This means that 7at_i is determined by Eq. (|l3| ) 
in (b) and (c) (therefore, in the vicinity of the outgoing 
resonance) and Eq. ([li]) is valid. Let us show that the 
r^r-i cancels out of the expression for the contribution 
of the diagram in Fig. [l| To do this note that Ajy^ from 
Eq. (|l5|) is proportional to the mean free path 
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(23) when £ N _ t » Xj, j = 0, 1, . . . , N - 2 . Thus, 
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where 

Tjv/3 is a dimetisionless function of Ao, Ai, . . . , Ajv_ 2 . For 
N = 2 we ixnS that (A^)" 1 = f°(z = 0) = 0, (Aai)" 1 = 
! x {z = 0) = 1/(A + A x ). This leads to T 20 = 0, T 21 = 1. 
Likewise, for N = 3, (A 3 oo) _1 = / 00 (z = 0) = 0, 
(Ash,)- 1 = I w (z = 0) = A /((A + Ai)(A + A 2 )), 
(Asn)- 1 = f n (z = 0) - Ai/((Ai + A )(Ai + A 2 )), 
(Asoi)- 1 - f 01 (z = 0) = A 2 /((A 2 + A )(A 2 + Ai)). For 
A 2 C 2 this lead to T 300 = 0, T 3W = A /(A + V), 
Tan = Xi/(Xo + Ai) and T 301 = 1. Using Eqs. ©, @ 
and (Eq) we obtain 



_! aujLom e 
2hK 



N-X 



(28) 



Thus, the quantity Tjv-i does not appear in the final 
result. 

Equations ^) and ( pl| ) for the contribution of Fig. ^ 
are valid when K z q is real which is true for all four fre- 
quency intervals. In (b), (c) and (d) the broadening 7^ _ 2 
is determined by the probability to emit an LO-phonon, 
whereas 7jv_ 2 — ► Tjv- 2 in (a), where Tjv- 2 << Jj, 
j = 0, 1, . . . , N - 3. Note that Eq. @) does not con- 
tent Tjv- 2 in the interval (a) as it was shown above. 
The contribution of Fig. || depends strongly on the be- 
havior of j eX c{K±, K z q) in the denominator of Eq. (|2l|). 
This is the inverse relaxation time of the exciton in 
the state with energy E exc (K±, K z o) — huji — (N — 
l)huj L o- For E exc (K^, K z0 ) > E 1H + huj LO the value 
of r Yexc(K±, Kzo) is determined by the probability to 
emit an LO-phonon and is proportional to a. How- 
ever, for E exc (K±, K z0 ) < E\h + h^LO, the real emis- 
sion of one LO-phonon is impossible and r y e xc(K±, K z q) 
is determined by other much weaker processes, so that 
Jexc(K±,Kzo) -> T exc (Kj_, K z0 ), with T exc « j exc . 
The change of the scattering mechanism dominating the 
broadening takes place at the frequency corresponding to 
the outgoing excitonic resonance. Below this point, the 
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contribution of Fig. || exceeds strongly the one of Fig. |[ 
Note that in this range we have to take into account other 
contributions involving processes with acoustic phonons 
(see below). 

Finally, the pole approximation (i.e., real transitions) 
for integrals over k$ z , k\ z , . . . , k^-\ z for the contribution 
of Fig. [l] and over ko z , ki z , . . . , fcjv-2z for Fig. || results 
in the constraints N > 2 and N > 3 for Eq. (^J) and 
Eq. (pi]), respectively. 



V. DISCUSSION AND CONCLUSIONS 

Let us consider the resonant behavior of the MPRRS 
efficiency as a function of H and loi. We limit ourselves 
to the case N > 3 where both Eq. (||) and Eq. ( pi] ) are 
valid. Both contributions increase in the vicinity of Kq = 
0, which corresponds to uJi m ax,m h ->oo — u g H + nu eH . 
Taking into account the finite value of leads to an 
exact relation LOi m ax{n) = lu 9 h +eHn/ /ic. This condition 
corresponds to the creation of EHPs in the vicinity of the 
Landau band bottoms. The resonant conditions can be 
achieved by changing either H or loi. The maxima in a 
magnetic field dependence take place at 



H„ 



\IC LOl — LO g 

7n+ 1/2 



(29) 



being independent on the order N of the scattering pro- 
cess. There is an additional resonanceQ for N = 3 corre- 
sponding to the contribution of Fig. || at 



+ neH/m e c + u>lo- This resonance follows from the 
increase of (A 2 i) _1 = 1/(A + Ai) = l/(hK /m e ^ + 
hKi/m e ji) in Eq. (^l|), when K\ — > 0, because of the di- 



vergence in 70 (see Eq. (h 
of Fig. || is equal to zerc 



I)). Note also that contribution 
3 for ui m in,m h ^oo(n, N - 1) = 
U g H + neH/m e c + (N — 1)lo LO - 

Above the outgoing resonance the contributions of 
Fig. and Fig. | in the MPRRS efficiency are of the 
same order of magnitude. However, as it was mentioned 
before, below the resonance the contribution of Fig. |^ 
strongly increases because of the strong increase in the 
exciton life time in the real intermediate state with the 
energy being too small for emission of an LO-phonon. In 
this range, other scattering processes like the absorption 
of LO-phonons, interaction with acoustic phonons, etc. 
have to be taken into account. We give now a qualita- 
tive picture of the process including into our considera- 
tion the distribution function of excitons with respect to 
K±, \K Z \. Let us introduce the integral efficiency for the 
N-th order process as 



Sn = 



^-dridcj, = — W s n 



dTLcLoe 



(30) 



where ui is the group velocity of incident light and 
W s n the normalized probability to emit the scattered 
light quantum per unit time.Q Equation (BOh differs from 



Eq. (||) only by the absence of the factor 4it5(loi — uo s — 
Nlolo)- 

On the other hand, 

X)WW= P e*cN-l( K ±A K *\M K ^\KS , (31) 

«s K±,K Z 

where P exc N-i(K±,\K z \) is the normalized dimension- 
less distribution function of excitons created by the light 
in a N-l LO-phonon-assisted process and ji(K±, \K Z \) 
the probability of an LO-phonon-assisted emission of 
the scattered light quantum which can be written as 
ji{K x , \K Z \) = V w s and 



2tt x - 



\p {f\Us\a){a\Hint\i) 
^ E t -E a + 



m-Ef) , (32) 



Hint is Frohlich interaction of the exciton with LO- 
phonons and U s represents the interaction of excitons 
with the light. The initial and final state energy is 
E{ = E exc (K±, \K Z \) and Ej = %lo s -{-Tilolo, respectively. 
The intermediate state energy E a = E a EHP + hu>LO in- 
cludes both the discrete and continuum part of the exci- 
tonic dispersion. Let us separate 7; in two corresponding 
parts, 7; = ~/idisc + licont- The outgoing resonance is re- 
lated to the contribution jidisco to jidisc coming from the 
transition via ground state of the exciton. According to 
Eq. @, we find 



7idisco{K±, \K Z \) = ( — ) |e s p c „| 2 

nc 5 \ m o / 



aoj 



l Q l exp (-a 2 H Kl/2) 



Kl + Kl ail a 2 H C 6 
[TijK^mh/M) - T){K zH m e lM)\ 2 

,(Kj_, l^.l) ~ UlH ~ ^LO? + llxcH (0)/4 



(33) 



Above the excitonic resonance, lo\ > lo\jj + Nlo^q, we 
have 

P IK WW W excN-l(Kj_, \K Z \) 

lexc{K^, \K Z \) 

where W exc N-i{K±,\K z \) is the normalized number of 
excitons created per unit time in the volume Vb in the 
N — 1 LO-phonon-assisted process. The probability 
W eX cN-i{K±,\K z \) has been calculated in Ref. H fo r 
N = 4. Being used in Eq. (UJ) together with Eqs. (|||), 
(3C) and ( |3l| ) it reproduces the result of Eqs. (||) and 

m 

Note that above the outgoing resonance the distrib 
tion is not zero only in very narrow interval of energie; 
since W exc N-i(K±, \K Z \) is proportional to 5[loi — (N — 
1)lolo — E exc (K±, \K Z \)/H]. However, for loi below the 
resonance the distribution becomes smooth. If the most 
important mechanism in this range is the interaction with 
acoustic phonons, one has to take into account diagrams 
with external acoustic phonon lines. In the range (b) 
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(see Eq. (Pq)) the smoothing of the distribution should 
be weaker than in the range (a). The reason of this is 
the kinetic energy of exciton in range (b) which is larger 
than the exciton binding energy. Since the probability of 
scattering and decay via the interaction with phonons are 
of the same order, the exciton decays after a few interac- 
tions with acoustic phonons. The decay of an exciton in 
the range (a) is suppressed because of its small energy. 
In this case, the distribution depends on the probability 
of non-radiative recombination. At zero magnetic field, 
the distribution of excitons in the range (a) has been 
considered in Refs. p|-^0| and p"T| . 

The smoothness of the distribution leads to (i) the 
broadening of the MPRRS peaks in the range (b) and 
especially in the range (a) and to (ii) the increase of 
the integral scattering intensity, since the diagrams with 
acoustic phonon lines give additional contributions into 
the MPRRS efficiency. 

To summarize, we have shown that the outgoing ex- 
citonic resonance has to be strongly enhanced under a 
high magnetic field. Above the outgoing resonance, the 
scattering efficiency for N > 3 may be up to a~ 2 times 
stronger than in a zero magnetic field where the MPRRS 
efficiencyu is proportional to a 3 , whereas in a high mag- 
netic field it is proportional to a, as it follows from 
Eqs. (||), (|^) and (HJ). The crossover from a 3 to a results 
from the quasi-one-dimcnsional character of free EHPs in 
N (Fig. |l|) or TV — 1 (Fig. ||) intermediate states under a 
high magnetic field. The enhancement is also valid for the 
ranges (a) and (b) below the excitonic resonance, where 
one has to calculate the exciton distribution function tak- 
ing into account the interaction with acoustic phonons. 
To the best of our knowledge such calculations have yet 
to be performed. However, the distribution function is 
proportional to the creation probability of excitons with 
energy in the interval between and Hlolo which is in- 
creased by a~ 2 times in a high magnetic field.d Thus, 
the MPRRS efficiency also increases below the excitonic 
resonance. 

The integral efficiency as a function of uii has to be 
asymmetric with respect to the point tkui = NTilulo + 
E\h because of the strong increase in the exciton life 
time below the resonance and appearance of additional 
contributions from the processes with acoustic phonons. 
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